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Apparatus and in vivo meth- 
ods for distinguishing normal and 
abnormal cervical tissue and de- 
tecting cervical intraepithelial neo- 
plasia (ON) in a diagnostic tissue 
sample. A source of electromag- 
netic radiation (36) is applied to 
the diagnostic cervical tissue sam- 
ple and the presumptively known 
normal cervical tissue. Fluores- 
cence intensity spectra induced in 
the diagnostic tissue sample and 
the presumptively known tissue by 
the radiation source is detected by 
at least one radiation detector (26). 
A computer (20) is programmed 
to measure a peak normal fluo- 
rescence intensity value in each 
f the normal fluorescence inten- 
sity spectra, to calculate an aver- 
age peak normal fluorescence in- 
tensity from the peak normal fluo- 
rescence intensity values, to mea- 
sure a peak fluorescence intensity 
value from the fluorescence inten- 
sity spectrum induced in the diag- 
nostic tissue sample, to calculate 
a slope parameter which is indica- 
tive of tissue abnormality, to cal- 
culate relative peak fluorescence . 
intensity as a function of the peak intensity value, the average peak normal intensity, and to detect tissue abnormality as a function 
f the slope parameter and the relative peak fluorescence intensity 
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DESCRIPTION 
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DIAGNOSIS OF DYSPLASIA USING 
LASER INDUCED FLUORESCENCE 

BACKGROUND OF THE INVENTION 

Field of the Invention 



The invention relates to methods and apparatus to 
differentiate between histologically normal and 
histologically abnormal tissues, and to differentiate 
neoplastic tissue from histologically abnormal non- 
15 neoplastic tissue. 

<;r:]-eenina for Cervical Intraepithelial Neo plasia fCINT 

Although there has been a significant decline in the 
20 incidence and mortality of invasive cervical carcinoma 
over the last 50 years, there has been an increase in 
both the reported and actual incidence of CIN. As a 
result, it has been estimated that the mortality of 
cervical carcinoma may rise by 20% in the years 2000-2004 
25 unless screening techniques for CIN are improved. 

Present screening for CIN and cervical cancer is 
relatively inexpensive but labor intensive because it 
initially relies on the results of a Pap smear; a false 

30 negative error rate of 20-30% is associated with 

insufficient cell sampling and/or inexpert reading of Pap 
smears. Given an abnormal Pap smear, colposcopic 
examination of the cervix (with a magnifying lens) 
followed by colposcopically directed biopsy and 

35 histologic examination of the tissue sample can provide a 
diagnosis of CIN. Histologic confirmation of the 
diagnosis, while relatively time-consuming and expensive, 
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is necessary because the accuracy of classification among 
abnormal tissues by colposcopy alone is limited, even in 
experienced hands. 

5 Improving the predictive value of colposcopy in 

distinguishing CIN from other abnormal tissues (e.g., 
tissue infected with human papilloma virus. (HPV) or 
inflammatory tissue) could reduce the number of required 
number of biopsies and thereby increase the speed and 

10 efficiency of the screening process. Diagnosis and 
treatment might be combined in a single office visit, 
with colposcopically directed treatments including loop 
electrosurgical procedures (LEEP) , cryo and laser 
therapies, and chemopreventive agents. Further, 

15 explication of improved methods to classify tissue as 
normal or abnormal, while not leading directly to a 
diagnosis of CIN, would reduce costs by allowing 
performance of colposcopy by medical technicians less 
skilled than trained colposcopists (usually physicians) . 

20 

Spectroscopic Methods i n Colposcopy 

Spectroscopic methods for differentiating cervical 
neoplasia from normal cervical tissue in vivo have been 

25 described. The methods rely generally on observations 

that the fluorescence of abnormal tissue is significantly 
weaker than that of normal tissue at several excitation 
wavelengths, e.g., 330, 350 and 450 nm. This property 
has been used for spectroscopic identification of 

30 histologically abnormal tissue. For example, in vitro 
fluorescence intensity comparisons at an excitation 
wavelength of 330 nm yielded positive predictive value, 
sensitivity and specificity of 86%, 88% and 75% 
respectively on colposcopically normal and abnormal 

35 biopsies from the same patient. Differences between 
neoplastic (CIN) and non-neoplastic abnormal tissues 
(inflammation and HPV infection) yielded the largest 
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spectrosc pic differences at an excitation wavelength of 
330 nm.. 

Additionally, fluorescence spectra have been 
5 measured in vivo to detect neoplastic tissues in 

different organ systems. A variety of methods for making 
such determinations have been proposed. For example, 
ratios of autof luorescence intensity at two different 
emission wavelengths have been used by many groups, and 
10 scores based on multi-variate linear or non-linear 

regressions and fits to extract concentrations of various 
chromophobes have been proposed by many others for 
inclusion in decision criteria. 

15 In one application, 337 nm wavelength excitation was 

applied to colonic tissue. Multi-variate linear 
regression analysis was used to correctly distinguish 
adenomatous polyps from normal colon and hyperplastic 
polyps with positive predictive value, sensitivity and 

20 specificity of 86%, 86% and 80% respectively. 

Tissue Classification Methods 

Previous attempts to reliably distinguish CIN from 
25 inflammation or HPV infection in vivo using colposcopy 

alone have been unsuccessful. Fluorescence intensity may 
be useful in this regard, but is not sufficient in itself 
because analogous fluorescence intensity measurements of 
the same tissue type may vary by more than a factor of 
30 two from patient to patient, and by about 15% within the 
same patient. Methods relating fluorescence intensities 
from abnormal and normal tissues of the same patient, 
however, tend to be more predictable and therefore more 
diagnostically useful. 

35 

Thus, in general, each patient must serve as h r own 
c ntrol. Considering cervical tissues in a giv n patient 
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excited with 337 nm wavelength electromagnetic 
excitation , tissues with HPV infection are less 
fluorescent with than tissues with chronic inflammation, 
and tissues with dysplastic changes exhibit even lower 
5 fluorescence than those with HPV infection. The lowest 
level of fluorescence (relative to analogous fluorescence 
measurements on other tissue types in the same patient) 
is exhibited by tissues with the most abnormal form of 
CIN. Analogous relationships among tissue fluorescence 
10 intensity measurements in a patient may exist if 

excitation wavelengths other than 337 nm are used because 
the shape and intensity of cervical tissue spectra do not 
change substantially when the excitation wavelength is 
increased or decreased by less than 10 nm. 

15 

Additional information useful for tissue 
classification may be found in the peak emission 
wavelength of tissues with CIN, which is positively 
correlated with the peak emission wavelengths of normal 
20 tissue spectra from the same patient. This relationship, 
however, is not observed for tissue samples with 
inflammation or HPV infection. 

Thus, a reliable method to spectroscopically 
25 classify tissue as normal or abnormal, and in the latter 
case to distinguish inflammation or HPV infection from 
CIN, is needed. Such a method could rely on one or more 
of the relationships described above, augmented with 
additional information indicative of the particular 
30 separation or classification desired. 

SUMMARY OF THE INVENTION 

35 The present invention includes in vivo spectroscopic 

methods and apparatus for differentiating normal tissue 
from abnormal tissue and for diagnosing CIN in diagnostic 
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cervical tissue samples. Diagnostic cervical tissue 
samples, are the tissue samples to b evaluated by the 
n n- invasive methods f th pres nt invention. Th 
methods include analysis of tissue fluorescence intensity 
5 spectra or portions thereof from both histologically 
normal tissue and histologically abnormal tissue in a 
patient desiring diagnosis. 

Reference in this application to histologically 
10 normal cervical tissue samples in a patient subject to 

the diagnostic methods of the present invention refers to 
tissue which is presumptively histologically normal. 
Normal tissue samples in such subjects are selected in 
vivo and tested non-invasively as described herein to 
15 ensure a substantial likelihood that they actually 
represent tissue which, if histologically evaluated, 
would be classified normal. 

In certain preferred embodiments, the spectra are 
20 represented on two-dimensional plots with fluorescence 
intensity being represented on the vertical axis and 
wavelength on the horizontal axis. For calculations 
involving slope measurements at predetermined wavelengths 
within spectra having different (unnormalized) peak 
25 values, each spectrum is normalized to its own maximum 

intensity value. For calculations involving the ratio of 
the peak intensity of a spectrum from unknown tissue 
(i.e., from a diagnostic cervical tissue sample) to the 
peak intensity of a spectrum from (presumptively) normal 
30 tissue in the same patient, normalization is not 
performed. 

Detecting Tissue Abnormality 

35 According the present invention, an in vivo method 

f detecting tissue abnormality in a diagnostic cervical 
tissue sample in a patient having known normal cervical 
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tissue comprises illuminating the diagnostic tissue 
sample and normal cervical tissue with electromagnetic 
radiation (preferably about 337 nm wavelength) . A 
plurality of normal fluorescence intensity spectra are 
5 detected from the known normal cervical tissue, the 
spectra being obtained serially or simultaneously, 
preferably through fiber optics which may be placed at a 
fixed distance from the cervical tissue. In some 
embodiments, spectra may be detected through analysis of 
10 a cervical image which includes the tissue areas to be 
sampled. 

A peak normal fluorescence intensity value is 
measured in each said normal fluorescence intensity 

15 spectrum, and an average peak normal fluorescence 

intensity calculated from said peak normal fluorescence 
intensity values. A fluorescence intensity spectrum is 
also detected from the diagnostic tissue sample, and a 
peak fluorescence intensity value measured from said 

20 fluorescence intensity spectrum. 

A predetermined portion of the fluorescence 
intensity spectrum (preferably in the region 
corresponding to wavelengths between about 410 nm and 430 

25 nm) from the diagnostic tissue sample furnishes slope 

information for calculation of a slope parameter which is 
indicative of tissue abnormality. The predetermined 
portion is empirically identified in prior clinical 
trials as facilitating accurate classification of tissue, 

3 0 and the slope parameter is preferably a function of 
average slope in the predetermined portion. Each 
spectrum is preferably normalized to its own maximum 
value prior to calculation of the slope parameter. 

35 Relative peak fluorescence intensity (preferably 

relative to peak normal fluorescence intensity) is 
calculated as a function of said peak fluorescence 
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intensity value, and tissue abnormality is detected as 
function of the calculated slope parameter and relative 
peak fluorescence intensity. The later function is 
preferably a predetermined empirical discriminant 
5 function obtained from prior clinical trials. The 

discriminant function may be either liner or nonlinear. 



n»». anting CIN 

10 According to the present invention, an in vivo 

method of detecting CIN in a diagnostic cervical tissue 
sample in a patient having known normal cervical tissue 
comprises classifying the diagnostic cervical tissue 
sample as abnormal by a method described herein, followed 

15 by illumination of the diagnostic tissue sample and known 
normal cervical tissue with electromagnetic radiation, 
preferably of about 337 nm wavelength. 



20 



A first fluorescence intensity spectrum from the 
diagnostic tissue sample is detected, as is a second . 
fluorescence intensity spectrum from the known normal 
cervical tissue. A first slope parameter which is 
indicative of CIN in the diagnostic cervical tissue 
sample is calculated from a predetermined portion of said 
25 first fluorescence intensity spectrum (corresponding to 
wavelengths between about 440 nm and 460 nm) . A second 
slope parameter which characterizes normal cervical 
tissue in the patient is calculated from a predetermined 
portion of said second fluorescence intensity spectrum 
30 (corresponding to wavelengths between about 410 nm and 
about 430 nm). Both first and second slope parameters 
are preferably calculated after normalization of the 
respective fluorescence intensity spectra to peak 
fluorescence intensity values of 1. Additionally, both 
35 first and second slope parameters are preferably 

functions of average slopes in said normalized first and 
sec nd fluorescence intensity spectra resp ctively in 
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regions lying substantially b tw en wavelengths 
corresponding to the respective predetermined portions. 

CIN is detected in the diagnostic cervical tissue 
sample as a function of said first and second slope 
parameters, the function preferably comprising a 
predetermined empirical discriminant function which, in 
some preferred embodiments is linear. 



BRIEF DESCRIPTION OF THE DRAWINGS 



Figure 1 is a schematic representation of apparatus to 
classify cervical tissue according to the present 
15 invention. 

Figure 2 is a schematic representation of apparatus for 
multi-pixel spatial-spectral imaging of the cervix. 

20 Figure 3 is a schematic representation of apparatus to 
acquire and process spectroscopic images of the 
cervix. 

Figure 4 represents a population of induced fluorescence 
25 intensity spectra representing colposcopically 

normal and histologically abnormal tissue, with a 
discriminant function for detecting tissue 
abnormality. 

3 0 Figure 5 represents a population of induced fluorescence 
intensity spectra representing histologically 
abnormal tissue (i.e., CIN, inflammation, HPV) , with 
a discriminant function for detecting CIN. 

35 Figures 6A-6F illustrate a flow chart to practice 

preferred embodiments of the methods of the present 
invention. 
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DETAILED DESCRIPTION 

CIN - cervical intraepithelial neoplasia 
5 HPV - human papilloma virus 

Normal /Abnormal Tissue Classification 

Normal/ abnormal tissue classification requires 
10 induced fluorescence intensity spectra from tissue areas 
known to be normal with a sufficiently high probability 
of accurate classification. Normal and abnormal cervical 
tissues are primarily identified in vivo by 
colposcopists. In any given patient, it is preferable to 
15 collect a plurality of fluorescence intensity spectra 

from tissue areas identified as normal by a colposcopist . 
An expected error rate of about 10-20% has been observed, 
which primarily represents samples which are 
colposcopically normal but in fact comprise inflammatory 
20 tissue. As explained below, provision is made for 

removal of erroneously- identified normal samples, with 
the result that those remaining are presumptively normal. 

As a patient is examined, peak normal fluorescence 
25 intensity values are measured on each spectrum from 

normal tissue. The intensity values are averaged and a 
standard deviation of normal peak values calculated. 
Depending on the degree of confidence required in the 
screening results and the skill of the operators, a 
30 portion of all spectra initially identified as 

colposcopically normal are discarded (i.e., the spectra 
are not included in calculation of the average peak 
normal fluorescence intensity, but may be assigned to the 
group of spectra to b classifi d subsequently as normal 
35 or abnormal) . Criteria f r discarding colposcopically 
normal samples may pr ferably require discarding all 
spectra associated with peak values which fall more than 
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one standard deviation below the average peak normal 
fluorescenc intensity. 

The peak fluorescence intensity of any test 
5 fluorescence spectrum (associated with cervical tissue to 
be evaluated for normality/abnormality) is divided by the 
corresponding average normal fluorescence intensity for 
that patient to yield a relative peak fluorescence 
intensity for that spectrum. Because normal tissue tends 

10 to have relatively high and uniform peak fluorescence 

intensities in any given patient, the above division will 
generally result in relative peak fluorescent intensities 
clustered around a value of approximately 1 for normal 
tissues. On the other hand, because abnormal tissue 

15 (whether characterized as inflammation, HPV infection, or 
CIN) tends to have lower-than-normal peak fluorescence 
intensities, the above division will generally result in 
relative peak fluorescent intensities clustered around a 
value substantially less than 1 for abnormal tissues. 

2 0 This condition furnishes a partial basis for classifying 
spectra as representing either normal or abnormal tissue. 

A" more complete basis for identification of spectra 
as representing either normal or abnormal tissue is 

25 provided by examination of a slope parameter associated 
with each spectrum, the parameter preferably being 
derived from slopes measured in the range 400-440 nm, 
preferably 410-430 nm, and in some preferred embodiments 
415-425 nm, as measured on a spectrum normalized to its 

30 own peak fluorescence value. The measured slope value is 
largely governed by two factors, i,e., the peak emission 
wavelength of fluorophoires contributing to the spectrum, 
and the reabsorption effect of oxy-hemoglobin with an 
absorption peak at 420 nm. As more hemoglobin 

35 contributes, the peak shifts to longer wavelengths and 
the slope increases. Similarly, as more NADH 
fluorescence contributes, the peak shifts to longer 
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wav 1 ngths and the slope increases. These effects are 
observed in diseased tissue, and probably account in part 
for the capacity f methods of the present invention to 
differentiate types of diseased tissue on the basis of 
5 their inducted fluorescence spectra. 



The parameter may be an average or tangential slope 
or other representative value of the range of actual 
slope values which will in general be found within the 

10 specified range (predetermined portion) in any normalized 
spectrum. Wavelengths defining the extent and location 
of the predetermined portion of a normalized spectrum 
which is indicative of tissue abnormality are determined 
from clinical trials involving fluorescence spectra 

15 obtained from tissue with histologically proven 

diagnoses. That is, the definition of the predetermined 
portion is empirically derived by comparing slopes over 
the entire wavelength range to find the portion of that 
range yielding the best discrimination using the methods 

20 described above. 

To evaluate a diagnostic cervical tissue sample, a 
relative peak fluorescence intensity and a slope 
parameter are calculated as explained above from its 

25 induced fluorescence intensity spectrum. The calculated 
values may then be plotted as a point on the graph of 
Figure 4 . The line 72 represents a predetermined 
empirical discriminant function which substantially 
separates points characterizing colposcopically normal 

30 and histologically abnormal tissue. If the point falls 
below line 72, the tissue sample is diagnosed as normal, 
whereas if the point falls above line 72, the tissue 
sample is diagnosed as abnormal. These decision criteria 
may be implemented mathematically as well as graphically 

35 by considering the equation f the lin which can be 

derived from Figure 4. Note that the position of line 72 
(slope and intercept) are subject to change with the 



WO 94/26168 



PCT/US94/05230 



-12- 

addition of more data points from clinical trials to 
those already pr sent on Figure 4. Note also, that the 
position of lin 72 may also be expected to change when 
populations having a substantially different incidence of 
5 abnormal cervical tissue from that illustrated in 
Figure 4 are considered. 

Differentiating ciN From Oth er Abnormal Tissues 

!0 Although normal/abnormal tissue classification on 

the basis of induced fluorescence spectra is useful for 
preliminary screening, patients identified as having 
abnormal cervical tissue must be further evaluated. In 
particular, CIN should be differentiated from 

15 inflammation or HPV infection because of the potential 
for CIN to progress to invasive cancer. 

As an aid to classification, a correlation has been 
observed between the wavelengths associated with peak 

20 intensity in corresponding fluorescence intensity spectra 
obtained from normal and neoplastic tissue in the same 
patient. Whereas three sample peak intensity wavelengths 
of normal spectra in a single patient were observed to be 
398, 426 or 442 nm, the corresponding peak intensity 

25 wavelengths from CIN spectra in the same patient were 

observed to be 442, 450 or 460, respectively. Thus the 
peak intensity of a fluorescence intensity spectrum tends 
to occur at longer wavelengths in fluorescence intensity 
spectra from tissue with CIN, compared to corresponding 

30 spectra from normal tissue in the same patient. This 
wavelength shift is relatively uncommon in abnormal 
tissues representing HPV infection or inflammation. 

A two-dimensional scatter plot was developed to map 
35 out the relationship in a population of patients b twe n 
the peak emission wavelength of th abnormal spectrum and 
average normal spectra from th same patient. This is 
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shown in Figure 5, where the abscissa corresponds to the 
slope of the abn rmal sp ctrum over the wavelength range 
440 nm to 460 nm, and the ordinate represents the average 
slope of the normal spectra from the same patient over 
5 the wavelength range 410-430 nm. All spectra were 
normalized to a peak intensity of 1 prior to slope 
calculation. Note that the slope of the spectra from 
samples with HPV infection or inflammation does not 
appear to correlate to the average slope of the 

10 corresponding normal spectra (linear correlation 

coefficient = 0.072). However, the slope of the spectra 
from samples with CIN displays a positive correlation 
(linear correlation coefficient - 0.442) to the average 
slope of the corresponding normal spectra. This enables 

15 differentiation of abnormal samples with CIN from those 
with HPV infection or inflammation. A predetermined 
empirical discriminant function which allows one to carry 
out the differentiation for each new patient is 
represented by the line 82 in Figure 5, which minimizes 

20 the number of misclassif ied samples. 

Thus, CIN tissue spectra are distinguishable from 
HPV or inflammation tissue spectra on a two-dimensional 
plot. The horizontal axis of the plot represents the 
25 value of a slope parameter obtained from an abnormal 

spectrum in the wavelength range of about 440 to 460 nm, 
while the vertical axis represents the value of a slope 
parameter obtained from a normal spectrum in the same 
patient in the wavelength range of about 410 to 430 nm. 

30 

A predetermined empirical function in 
two-dimensional space substantially separates points 
corresponding to CIN tissue fluorescence spectra from 
points corresponding to HPV or inflammation tissue 
35 sp ctra. In c rtain preferred embodiments, this function 
may be linear and acceptably minimize the number of 
misclassif ied points. In other preferred embodiments, 
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the function may describe a nonlinear decision surface. 
Determination of the preferred decision surface for any 
population depends on the degree to which the population 
is characterized by clinical data used to estimate the 
5 decision surface. 

Apparatus for tissue clas sification 

Figure 1 is a schematic representation of apparatus 

10 to classify cervical tissue according to the present 
invention. To obtain tissue fluorescence spectra, 
electromagnetic radiation (e.g., light, in certain 
embodiments) in the form of laser light from nitrogen 
laser 36 is applied through coupling lens 34 and single 

15 fiber optic excitation fiber 32 to probe 38. Also within 
probe 38 are collection fiber optic fibers 30 and quartz 
shield 40 , the shield 40 acting to keep fibers 30 and 32 
properly spaced from any surface to which probe 40 is 
applied, application to a cervix preferably being under 

20 colposcopic observation. Resulting tissue fluorescence 

is transmitted by fibers 30 through coupling lenses 28 to 
polychromator 27 , and thence to intensified diode array 
26. Array 26 , controlled by controller 22 through gate 
pulser 24, detects fluorescence intensity spectra which 

25 are relayed through controller 22 to computer 20. 
Computer 20 is programmed to classify tissue in 
accordance with methods described herein. 

Note that laser light is not necessary for practice 
30 of the present invention. Laser 36 may be replaced in 
some embodiments by an incandescent or other type lamp 
with an associated f ilter to produce quasi-monochromatic 
light. Additionally, fibers 32 and 30 may be replaced in 
some embodiments with fibers serving both the function of 
35 illumination (excitation) and transmission of tissue 
fluorescence. Intensified diode array 26 and 
polychromator 27 may b r placed by a subassembly 
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comprising radiation filters and photomultiplier tubes to 
reduce costs in certain embodiments. 

Figure 2 schematically illustrates apparatus used in 
5 certain embodiments of the present invention to 

simultaneously collect fluorescence spectra from multiple 
areas of the cervix. A source 50 of electromagnetic 
radiation (delivered through fiber or non-fiber optics) 
is used to illuminate the desired areas of the cervix 

10 (sample 51) , including both some normal and some abnormal 
areas. A geometric array 53 of fibers 52 collects tissue 
'fluorescence originating from known normal and unknown 
regions of the cervix, the latter to be diagnosed. The 
ends of fibers 52 proximate to imaging polychromator 54 

15 are arranged in a linear array 55 at the entrance slit to 
polychromator 54. Polychromator 54 is coupled to a 
charge coupled device camera 56. Polychromator 54 
disperses wavelength across one axis of the array, and 
position on the cervix varies across the other dimension. 

20 Thus this system provides a spatial-spectral image on 
computer-display 58 of fluorescence spectroscopic 
information in the cervix. 

The operator identifies one or more fibers which 
25 view normal cervix. The fluorescence intensity spectrum 
from each fiber viewing an unknown area of the cervix is 
then processed by methods described herein to determine 
whether the tissue is histologically abnormal and whether 
CIN is present. This information can be presented as a 
30 spatial image of tissue histologic condition on computer- 
display 58. 

Figure 3 is a schematic illustration of apparatus 
for certain embodiments of the present invention. As in 
35 Figure 2, illumination of sample 51 with electromagnetic 
radiation from s urce 50 results in flu rescenc 
intensity spectra which are sensed by charge coupled 
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device 56 after passing through variabl band pass filter 
62 and colposcope imaging optics 60. Computer/ display 
unit 68 accepts spectra from charge coupled d vice 56 for 
processing and display. 

5 

This system provides a two-dimensional image of 
variation in fluorescence intensity as a function of 
position on the cervix at specific wavelength band 
governed by the transmission characteristics (optical 

10 properties) of imaging optics 60 and filter 62, which can 
be rapidly changed. Images are acquired at bands 
centered near 398, 410, 426, 442, 450 and 460 nm 
sequentially. The maximum intensity I MAX is identified 
for each pixel. A composite image indicating the peak 

15 intensity as a function of position is then formed from 
these images. Images indicating the slopes of the 
fluorescence spectra from the tissue in each pixel over 
the ranges 410-430 nm and 440-460 nm are calculated 
according to 1(426) - I (410) /[ (426-410) I MAX ] and 1(460) 

20 - 1(442) /[ (460-442) I MAX ]. The operator identifies one 
or more pixels corresponding to normal regions of the 
cervix. A relative intensity image is constructed by 
dividing the peak intensity image by the average normal 
intensity. Thus, for each pixel the relative intensity 

25 and slopes at 410-430 nm and 440-460 nm are available. 
These data are used to classify the state of the tissue 
in each pixel according to the methods presented herein. 
This information can be presented on computer-display 
unit 68 as a spatial image of tissue histologic 

3 0 condition. 

Obtaining Decision Surfaces 

Data points in Figure 4 represent induced 
35 fluorescence intensity spectra obtain d from cervical 
tissue in a population of patients, each patient having 
both colposcopically normal tissue and histologically 
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abnormal tissue. The latter tissue includes CIN, 
inflammatory tissue, and tissu having HPV infection. 
Note that fluorescence intensities plotted along the 
abscissa of Figure 4 are relative to (divided by) average 
5 peak normal fluorescence intensity obtained from the 

population of spectra considered. Spectra from the two 
tissue groups are processed in accordance with methods 
described herein to classify the respective tissues as 
normal or abnormal. Decision surface 72 (a line in 

10 2-space) is empirically established to minimize overlap 
of substantially normal 76 and substantially abnormal 74 
groups. The equation of surface 72 as determined from 
Figure 4 constitutes a predetermined discriminant 
function useful in detecting tissue abnormality in any 

15 diagnostic cervical tissue sample when fluorescence 

intensity spectra are induced in the tissue sample and 
processed according to methods described herein. 

Data points in Figure 5 represent induced 

20 fluorescence intensity spectra obtained from cervical 

tissue in a population containing histologically abnormal 
tissue, the latter tissue including CIN, inflammatory 
' tissue, and tissue having HPV infection. Spectra from 
the three tissue groups are processed in accordance with 

25 methods described herein to classify the respective 

tissues as characteristic of CIN or not characteristic of 
CIN. Decision surface 82 (a line in 2-space) is 
empirically established to minimize overlap of 
substantially CIN 86 and substantially not-CIN 84 groups. 

30 The equation of surface 82 as determined from Figure 5 
constitutes a predetermined discriminant function useful 
/in detecting CIN in any histologically abnormal 
diagnostic cervical tissue sample when fluorescence 
intensity spectra are induced in th tissue sample and 

35 processed according to the methods described herein. 
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Figur s 6A-6F illustrate a flow chart to practice 
preferred embodiments of the methods of the present 
invention. Note that steps in the flow chart continue 
from Figure 6A to Figure 6B f from Figure 6B to Figure 6C, 
5 and from Figure 6C to Figure 6D. Figures 6E and 6F 

represent subroutines to apply the particular decision 
rules applicable to either separation of tissue into 
normal and abnormal classifications (Figure 6E) , or 
identification of CIN in samples of abnormal tissue 

10 (which may include CIN, inflammation or HPV infection) 

(Figure 6F) . The subroutines (subcharts) illustrated in 
Figures 6E and 6F are called for within the steps 
illustrated in Figure 6D. It is assumed in applying the 
steps called for in Figures 6A-6F that sufficient spectra 

15 are made available to ensure that normal cervical tissue 
in the patient is identified and characterized by 
fluorescence spectra with a sufficiently high likelihood 
of accurate classification, as indicated in the methods 
described herein. 

20 Referring to Figs. 6A-6F, the process steps executed 

in the diagnostic method of the present invention are 
shown. After the method is started, the following 
parameters are initialized in step 101: 

25 counters i, 1c, I are initialized to 1; 

N Q - the number of normal spectra to acquire; 
N u - the number of unknown spectra to acquire; 
Mj - the slope of the normal /abnormal decision 
surface = 0.0088; 
3 0 B l " the intercept of the normal abnormal decision 

surface = -0.0005; 
M 2 - the slope of the CIN/not CIN decision surface = 

4.76; 

B 2 - the intercept of the CIN/not CIN decision 
35 surface = 0.0045. 

In step 102, the operator identifies the ith normal 
site on the cervix. The ith normal spectrum is acquired 
and stored in memory as N}(Xj) in step 103. In step 104, 
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the largest (peak) intensity valu is extract d from 
Nj(Xj); this value is stor d in memory as P { . In step 105, 
the normalized normal spectrum from site i (NNj(Xj)) is 
calculated as NNj(Xj)=Nj(Xj) /P| and stored in memory. In 
5 step 106, the slope of the spectrum from 410-430 nm, Sj, 
is calculated from NNj(Xj) as Sp[NNj(430)-NNj(410) ]/ (430 nm 
- 410 nm) . This value is stored in memory in step 106. 

In step 107 the value of i is compared to N Q . If i 
10 f N 0 , then i is reset to i+1 in step 108 and the method 
returns to step 102. Steps 102 through 107 are repeated 
until counter i is equal to N Q . At this time, step 109 
is executed, wherein the average peak intensity of normal 

spectra, P" , is calculated as P=V Pt/N 0 and stored in 

2=1 

memory. Step 110 is then executed in which the standard 
deviation of the peak intensity, a, is calculated as 

15 



0 = 



£ (P r P) L j (N 0 -l) . This value is stored in memory 



and then block 111 is executed wherein counters i, £ and 
k are reset to 1. 
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Blocks 112 through 122 determine which of the 
spectra collected from areas identified as normal will be 
considered normal within the method and which will later 
be reclassified as normal or abnormal in later steps of 
5 the method. In block 112, the values of Pj and P are 
compared. If - P is greater than a the spectrum from 
the ith sample (N:(X:)) is renamed and stored in memory as 
U k (Xj) in block 113. The type of this sample is 
classified later within the method. Following this, 

10 blocks 114 and 115 are executed where the counters k and 
i are incremented by one respectively. Control is then 
returned to block 112 of the method. If Pj ■- 1? is not 
greater than o, the peak intensity of the corresponding 
spectrum (Pj) is renamed and stored as Pj in block 116. 

15 In block 117, executed next, the slope of the spectrum 
from 410-430 nm (Sj) is renamed and stored as S^. The 
type of this sample is considered normal within the 
method. In block 118 the method tests to see whether i = 
N Q . If not, in block 119, the counter I is incremented 

2 0 by one and in block 115 the counter i is incremented by • 
one. Control is returned to block 112 and the sequence 
from 112 to 118 is repeated until all acquired spectra 
have been tested and i is equal to N Q . At this time, 
control is shifted to block 120 and N Q is reset to I and 

25 stored in memory. 
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In block 121, the average peak intensity of those 
samples considered normal is r calculated (P) and stored 

according to P= £ Pf/ty> In block 122 the average slope 

£=l 

of those samples considered normal (S) is calculated and 

5 stored in memory according to 5= £ S(/N 0 . 

£=l 



In block 123, the operator is directed to acquire 
the kth unknown spectru, U k (Xj) ; this is stored in memory. 
In block 124 it is determined whether the desired number 

10 of unknown spectra, N u , have been collected by testing 

whether k ■ N„. If not, in block 125 k is incremented by 
one and control is returned to block 123. The steps in 
blocks 123 to 124 are repeated until k = N u , at which 
time, control is shifted to block 126 and k is 

15 reinitialized to 1. 

Control is then passed to block 127, where the peak 
intensity is extracted from the kth unknown spectrum, 
U k (Xj), and stored as UP k . The relative peak intensity of 
20 the kth unknown spectrum, UR k , is then calculated as UR k = 
UP k /P* in block 128 and stored. The normalized unknown 
spectrum, NU k (Xj) is then calculated as NU k (Xj) =U k (Xj) /UP k 
and stor d in block 129. In block 130, the slope of the 
kth normalized unknown spectrum from 410-430 nm (US lk ) is 
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calculated as US lk =(NU k (430nm)-NU k (410 nm))/(430 nm-410 
ran) and stored. 

In block 132 the type of the kth unknown sample is 
5 determined as normal or abnormal by following the steps 
in Figure 6E. In block 132a, the parameter Y kl is 
calculated for the kth unknown spectrum as Y kl - MjURjj + 
Bj; this value is stored in memory. In block 132b, US lk 
is compared to Y kl . If US lk is less than Y kl , the sample 

10 is classified as normal; this information is stored in 
memory in block 133 and control is transferred to block 
137. If US lk is not less than*Y kl , the sample is 
classified as abnormal and control is passed to block 
134. In block 134, the steps of Figure 6F are used to 

15 diagnose whether or not the abnormal sample contains CIN. 
In block 134a, the parameter Yj^ is calculated for the kth 
unknown spectrum as Yj^ = M 2 US 2k +B 2 . This value is stored 
in memory. In block 134b, the values of S and Y^ are 
compared. If S is less than Y^ then the sample is 

20 classified as CIN and this information is stored in 

memory in block 135. If S is not less than Y^ then the 
sample is classified as abnormal, but not containing CIN, 
and this information is stored in memory in block 136. 
When the sample type of the kth unknown sample has been 

25 classified, in blocks 133, 135 or 136, control is passed 
to block 137, where the values of k and N u are compared. 
If k is not equal to N„, the counter k is incremented by 
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one in block 138 and control returns to block 127. The 
steps in bl cks 127 through 137 are repeated until k is 
equal to N u , indicating that all unknown samples have 
been classified and the method is stopped in block 139. 
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CLAIMS 



5 



1. An in vivo method of detecting tissue abnormality in 
a diagnostic cervical tissue sample, comprising: 

illuminating the diagnostic tissue sample with 
electromagnetic radiation; 



detecting a fluorescence intensity spectrum from the 
10 diagnostic tissue sample; 

measuring a peak diagnostic tissue fluorescence 
intensity value from said fluorescence 
intensity spectrum; 



15 



20 



calculating from a predetermined portion of said 
fluorescence intensity spectrum a diagnostic 
tissue slope parameter which is indicative of 
tissue abnormality; 

calculating relative peak diagnostic tissue 

fluorescence intensity as a function of said 
peak fluorescence intensity value; and 



25 detecting tissue abnormality as a function of said 

diagnostic tissue slope parameter and said 
relative peak diagnostic tissue fluorescence 
intensity. 
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2. The method of claim 1, said illuminating step 
comprising, illuminating said diagnostic tissue sample 
with electromagnetic radiation having a wavelength of 

5 about 337 nm. 

3. The method of claim 1, said step of calculating said 
diagnostic tissue slope parameter comprising,- calculating 

10 said diagnostic tissue slope parameter from a portion of 
said fluorescence intensity spectrum corresponding to 
wavelengths between about 410 nm and about 4 30 nm. 

15 4, The method of claim 1, wherein said slope parameter 
is calculated after normalization of said fluorescence 
intensity spectrum with said peak fluorescence intensity 
value. 

20 

5. The method of claim 4, said step of calculating said 
diagnostic tissue slope parameter comprising, calculating 
said diagnostic tissue slope parameter as a function of 
an average slope in said normalized fluorescence 

25 intensity spectrum in said predetermined portion. 

6. The method of claim 1, further comprising 
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illuminating a known normal cervical tissue with 
said electromagnetic radiation; 



detecting a plurality of normal fluorescence 
5 intensity spectra from the known normal 

cervical tissue; 

measuring a peak normal tissue fluorescence 
intensity value in each said normal 
10 fluorescence intensity spectrum; and 

calculating an average peak normal tissue 

fluorescence intensity from said peak normal 
tissue fluorescence intensity values; 

15 

said step of calculating relative peak fluorescence 
diagnostic tissue fluorescence intensity 
comprising, calculating relative peak 
diagnostic tissue fluorescence intensity as a 
20 function of both said peak fluorescence 

intensity value and said average peak normal 
tissue fluorescence intensity. 



25 7. The method of claim 6, wherein said calculating 
relative peak fluorescence intensity step comprises 
calculating relative peak fluorescence intensity as a 
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f unction of both said peak flu rescence int nsity value 
and said average peak normal fluorescence intensity. 



5 8. The method of claim 1, said step of detecting tissue 
abnormality comprising, detecting tissue abnormality as a 
predetermined empirical discriminant function of said 
abnormal tissue slope parameter and said relative peak 
abnormal tissue fluorescence intensity. 



9. The method of claim 8, said step of detecting tissue 
abnormality comprising, detecting tissue abnormality as a 
predetermined linear empirical discriminant function of 
15 said abnormal tissue slope parameter and said relative 
peak abnormal tissue fluorescence intensity. 



10. An in vivo method of detecting tissue abnormality in 
20 a diagnostic cervical tissue sample in a patient, 
comprising: 

identifying a presumptively histologically normal 
cervical tissue sample from the patient; 



25 



illuminating th diagnostic tissu sample and said 
presumptiv ly histologically normal tissue 
sample with lectromagn tic radiation; 



WO 94/261*8 PCT/US94/05230 

-28- 

detecting a first fluorescence intensity spectrum 
from the diagnostic tissu sample; 

detecting a second fluorescence intensity spectrum 
5 . from said presumptively histologically normal 

tissue sample; 

calculating a slope parameter indicative of tissue 
abnormality from said first fluorescence 
10 intensity spectrum; 

calculating an intensity parameter characteristic of 
normal tissue from said second fluorescence 
intensity spectrum; 

15 

calculating a discriminant parameter as a function 
of said slope and intensity parameters; 

detecting tissue abnormality in the diagnostic 
20 cervical tissue sample by substituting said 

discriminant parameter in a predetermined 
discriminant function. 



25 



11. An in vivo method of detecting CIN in a diagnostic 
cervical tissue sample in a patient having presumptively 
known normal cervical tissue, comprising: 
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illuminating the diagnostic tissue sample and 

presumptively kn wn normal cervical tissue with 
electromagnetic radiati n; 

5 detecting a first fluorescence intensity spectrum 

from the diagnostic tissue sample; 

detecting a second fluorescence intensity spectrum 
from the -presumptively known normal cervical 
10 tissue; 



calculating from a predetermined portion of said 

first fluorescence intensity spectrum a first 
slope parameter which is indicative of CIN in 
15 the diagnostic cervical tissue sample; 

calculating from a predetermined portion of said 

second fluorescence intensity spectrum a second 
slope parameter which characterizes 
20 presumptively normal cervical tissue in the 

patient; and 



25 



detecting CIN in the diagnostic cervical tissue 

sample as a function of said first and second 
slope parameters. 
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12. The method of claim 11, said illuminating step 
comprising, illuminating said diagnostic tissue sample 
with electromagnetic radiation having a wavelength of 
about 337 nm. 

5 

13. The method of claim 11, said step of calculating 
said first slope parameter comprising, calculating said 
first slope parameter from a portion of said first 

10 fluorescence intensity spectrum corresponding to 
wavelengths between about 440 nm and about 460 nm. 

14. The method of claim 11, said step of calculating 

15 said second slope parameter comprising, calculating said 
second slope parameter from a portion of said second 
fluorescence intensity spectrum corresponding to 
wavelengths between about 410 nm and about 430 nm. 

20 

15. The method of claim 11, wherein said first and 
second slope parameters are calculated after 
normalization of said first and second fluorescence 
intensity spectra to peak fluorescence intensity values 

25 of 1. 
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16. The m thod of claim 15, wherein said first and 
second slope parameters are functions of average slopes 
in said normalized first and second f luorescenc 
intensity spectra respectively, said slopes being 
5 -measured in predetermined regions of said normalized 
first and second fluorescence intensity spectra 
respectively, said regions lying substantially between 
wavelengths corresponding to said predetermined portions 
of said first and second fluorescence intensity spectra 
10 respectively. 



17. The method of claim 11, said step of detecting CIN 
comprising, detecting CIN as a predetermined empirical 

15 discriminant function of said first and second slope 
parameters . 

18. The method of claim 17, said step of detecting CIN 
20 comprising, detecting CIN as a predetermined linear 

empirical discriminant function of said first and second 
slope parameters. 



25 19. An apparatus for in vivo detection of tissue 

abnormality in a diagnostic cervical tissue sample in a 
patient having presumptively kn wn normal cervical 
tissue, comprising: 
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source of electromagnetic radiation to be applied 
to the diagnostic cervical tissue sample and 
the presumptively known normal cervical tissue; 

least one radiation detector adapted to detect a 
fluorescence intensity spectrum induced in the 
diagnostic cervical tissue sample and a 
plurality of normal fluorescence intensity 
spectra induced in the presumptively known 
cervical tissue by said radiation source; 

computer, connected to said source of 

electromagnetic radiation and said radiation 
detector, programmed to measure a peak normal 
fluorescence intensity value in each said 
normal fluorescence intensity spectrum, to 
calculate an average peak normal fluorescence 
intensity from said peak normal fluorescence 
intensity values , to measure a peak 
fluorescence intensity value from said 
fluorescence intensity spectrum induced in the 
diagnostic cervical tissue sample, to calculate 
from a predetermined portion of said 
fluorescence intensity spectrum a slope 
parameter which is indicative of tissue 
abnormality, to calculate relative peak 
fluorescence intensity as a function of said 
p ak fluorescence int nsity value and said 
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average peak normal fluorescence intensity, and 
t detect tissue abnormality as a function of 
said slope parameter and said relative peak 
fluorescence intensity. 



WO 94/26168 



1/11 




FIGURE 1 
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FIGURE 3 
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START: i = k = / = l 

READ- No = NUMBER OF NORMAL SPECTRA TO ACQUIRE 
N u = NUMBER OF UNKNOWN SPECTRA TO ACQUIRE 
Ml, Bi = SLOPE AND INTERCEPT OF NORMAL/ABNORMAL 

DECISION SURFACE 
M2, B2 = SLOPE AND INTERCEPT OF CIN/ NOT C1N 
DECISION SURFACE 
(Mi = 0.0088, Bi = -0.0005, M2 = 4.76, B2 = 0.0045) 

101 



T 



IDENTIFY NORMAL SITE ON CERVIX 102 1 



I 



ACQUIRE & NORMAL SPECTRUM; STORE AS Nj (Xj) 103 



I 



EXTRACT PEAK INTE NSITY; STORE AS Pj 

1 



104 



CALCULATE AND STORE NORMALIZED NORMAL SPECTRUM; 
NNi(Xj)=Ni(Xj)/Pi 105 



EXTRACT SLOPES FROM (410-430) NM FROM 

NNj(X,); STORE AS Sj 106 




NO 



10* 
i = i+l - 



YES 



© 



FIGURE 6A 
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CALCULATE AVERAGE PEAK INTENSITY: 
N 0 

P = SPi/N 0 



i=l 



109 



CALCULATE STANDARD DEVIATION: 
N 0 

ff = SQRT(S(P i -P) 2 /(N 0 -l)) 
i=l 



110 



, = 1, 1=1, k=l 



111 



i = i + l 




STORE SPECTRUM AS Ufc: 
Uk(Xj)=Ni(X.j) 



k=k+l - 



STORE PEAK INTENSITY: P/ = P 



116 




FTOTTRF. 6R 
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RE-CALCULATE AVERAGE PEAK INTENSITY: 
N 0 

F = Epf/N 0 



121 



CALCULATE AVERAGE SLOPE: 

S = Ss//N 0 
l = i 



122 



ACQUIRE ^ UNKNOWN SPECTRUM; 
STORE AS Ufc(Xj) 




FIGURE 6C 
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EXTRACT PEAK INTENSITY FROM U k 
STORE AS UPfc 



CALCULATE RELATIVE PEAK INTENSITY: 

UR k = UP k / P 12g 



CALCULATE AND STORE NORMALIZED 
UNKNOWN SPECTRUM: 

NU k (Xj)=U k (XjyUP k 129 



EXTRACT SLOPES FROM (410-430) NM FROM 
NU k (Xj); STORE AS US ik 



EXTRACT SLOPES FROM (440-460) NM FROM 
NU k (Aj); STORE AS US2k 



DIAGNOSE AS NORMAL OR ABNORMAL 
(SEE SUB-CHART 6E) 



132 



133- 



138 



NORMAL 
"(STORE) 



137- 



ABNORMAL 



DIAGNOSE AS 
CIN/NOTCIN 
(SEE SUB-CHART 6F) 



134 



CIN 
(STORE) 



NOT CIN 
(STORE) 





k = k + l 


NO 


4 






136 



FIGURE 6D 
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CALCULATE Ykl = Mi UR k + B] ^ 






ISUSi k <Y k i n2b 




— ►ABNORMAL 




YES 


1 


r 



NORMAL 



FIGURE 6E 
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CALCULATE Y& = M 2 US 2 k + B 2 134a 






IS S<Y k2 134b 


a 


— ► NOT ON 

YES 

r 



FIGURE 6F 
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